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ABSTRACT 



We report the detection of the cool, Jovian-mass planet MOA-2007-BLG- 
400Lb. The planet was detected in a high-magnification microlensing event (with 
peak magnification A max = 628) in which the primary lens transited the source, 
resulting in a dramatic smoothing of the peak of the event. The angular extent 
of the region of perturbation due to the planet is significantly smaller than the 
angular size of the source, and as a result the planetary signature is also smoothed 
out by the finite source size. Thus the deviation from a single-lens fit is broad and 
relatively weak (~ few percent). Nevertheless, we demonstrate that the planetary 
nature of the deviation can be unambiguously ascertained from the gross features 
of the residuals, and detailed analysis yields a fairly precise planet /star mass ratio 
of q = 2.6 ± 0.4 x 10~ 3 , in accord with the large significance (A% 2 = 1070) of the 
detection. The planet/star projected separation is subject to a strong close/wide 
degeneracy, leading to two indistinguishable solutions that differ in separation 
by a factor of ~ 8.5. Upper limits on flux from the lens constrain its mass to be 
M < 0.75 M Q (assuming it is a main-sequence star). A Bayesian analysis that 
includes all available observational constraints indicates a primary in the Galactic 
bulge with a mass of ~ 0.2 — O.5M and thus a planet mass of ~ 0.5 — 1.3M Jup . 
The separation and equilibrium temperature are ~ 0.6 — 1.1 AU (~ 5.3 — 9.7 AU) 
and ~ 103 K (~ 34 K) for the close (wide) solution. If the primary is a main- 
sequence star, follow-up observations would enable the detection of its light and 
so a measurement of its mass and distance. 



1. Introduction 

In the currently favored paradigm of planet formation, the location of the snow line in 
the protoplanetary disk plays a pivotal role. Beyond the snow line, ices can condense, and 
the surface density of solids is expected to be higher by a factor of several relative to its 
value just inside this line. As a result of this increased surface density, planet formation is 
expected to be most efficient just beyond the snow line, whereas for increasing distances from 
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the central star the planet for mation efficien cy drops, as the surface density decreases and 
the dynamical time increases (jLissauerl 119871 ). In this scenario, gas-giant planets must form 
in the region of the protoplanetary disk immediately beyond the snow-line, as the higher 
surface density is required to build icy protoplanetary cores that are sufficien tly massive to 



accre te a substantial gaseous envelope while there is remaining nebular gas (jPollack et al. 



19961 ). Low-mass primaries are expected to be much less efficient at forming gas giants 
becau se of the longer dynamical times and lower surface densities at th e snow lines of these 
stars (ILaughlin et al.l 120041 ; llda fc Lml 120051 ; Kennedy fc Kenyonl 120081 ) . Migration due to 
nebular tides and other dynamical processes can then bring the icy c ores or gas giants from 
their formation sites to o rbits substantially interior to the snow line (ILin et al.lll996l ; IWard 
1997l : lRasio fc Fordlll996h . 



The precise locatio n of the snow line in protoplanetary disks is a matter of some debate 



^e.g., iLecar et al.ll2006l ). and is even likely to evolve during the epoch of planet formation, 
particularly for low- mass stars ([Kennedy et al.ll2006l ; Kennedy fc Kenyonl 120081 ) . The con- 
densation temperature of water is ~ 170 K, and a fiducial value for the location of the snow 
line in solar-mass stars motivated by observations in our solar system is ~ 2.7 AU. This may 
scale linearly with the stellar mass M, since the ste llar luminosity during the e poch of planet 
formation scales as ~ M 2 for stars with M < M (IBurrows et al.lll993l . 119971 ). Whereas the 
radial velocity and especially transit methods are most sensitive to planets that are close to 
their parent star at distances well inside the snow line, the sensitivity of the microlensing 
method peaks at planetary separations near the Einstein ring radius of the primary lens 
jMao fc Paczvnskilll99ll : bould fc Loe~blll992h . which is ~ 3.5 AU(M/M ) 1 / 2 for typical lens 
and source distances of 6 kpc and 8 kpc, respectively. This corresponds to a peak sensitivity 
at equilibrium temperatures of T eq ~ 150 K(M/M Q ) for a mass-luminosity relation of the 
form L oc M 5 , and distances relative to the snow line of ~ 1.3(M/M Q )~ 1/2 if the location of 
the snow line at the epoch of planet formation scales as M. Thus microlensing is currently 
the best method of prob ing planetary systems in the critical region just beyond the snow 
line f lGould fc Loeblll992h . 



Planetary perturbations in microlensing events come in two general classes. The major- 
ity of planetary perturbations are expected to occur when a planet directly perturbs one of 
the two images cre ated by the primary lens, as the image sweeps by the planet during the 
microlensing event (IGould &: Loeblll992l ). Although these perturbations are more common, 
they are also unpredictable and can occur at any time during the event. Early microlensing 
planet searches focused on this class of perturbations, as it was the first to be identified an d 
explored theoretically JGould fc Loebl Il992l ; bennett & Rhiel Il996l : baudi & Gouldl Il997h . 
The second class of planetary perturbations occurs in high- magnification events, in w hich 
the source becomes very closely aligned with the primary lens ( IGriest fc Safizadehlll998l ). In 
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such events, the two primary- lens images b ecome highly distorted and sweep along nearly the 
entirety of the Einstein ring (jLiebea Il964l ) . These sweeping images probe subtle distortions 
of the Einstein ring caused by nearby p lanets, which will give rise to perturbation s within the 
full- width half-maximum of the event ( Bond et al. 2002 ; Rattenbury et al]|2002l ). Although 
high-magnification events are rare and so contribute a minority of the planetary perturba- 
tions, they are individually more sensitive to planets because the images probe nearly the 
entire Einstein ring. Furthermore, since the perturbations are localized to the peak of the 
event which can be predicted beforehand, they can be monitored more efficiently with limited 
resources than the more common low-magnification events. 

For these reasons, current microlensing planet searches tend to deliberately focus on 
high-magnification events. Thus , of the seven prior microlensing planets discovered to date 
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20081 ). five have been found in high- magnification events, with peak mag- 
nifications ranging from A = 40 to A = 800. However, despite the fact that this planet-search 
strategy has proven to be so successful, the properties of the planetary perturbations gener- 
ated in high-magnification events are less well-understood than those in low-magnification 
events. 

Most of the studies of the properties of planetary perturbations in high-magnification 
events have focused on the properties of the caustics, the locus of points defining one or more 
closed curves, upon which the magnification of a point source is formally infinite. The mor- 
phology and extent of the region of significant perturbation by the planetary companion can 
be largely understood by the shape and size of these caustic curves. Planetary perturbations 
in high-magnification events are caused by a central caustic located near the position of the 
primary, and thus several authors have considered the size and shape of these central caus 



tics as a function of the parameters of the planet (IGriest fc Safizadehl Il998l ; iDominiki 11999 



Chung et al.ll2005l ). These and other authors have identified several potential degeneracies 
that complicate the unique interpretation of central caustic perturbations. The first to be 
identified is a degeneracy such that the caustic structure (and so light curve morphology) of a 
planet with mass ratio q <C 1 and projected separation in units of the Ei nstein ring d not too 
close to unity is essentially identical under the transformation d <-» d^ 1 (IGriest fc Safizadeh 
19981 ). A second degeneracy arises from the fact that very close or very wide roughly equal- 
mass binaries also produce perturbations near the peak of the light curves. These give rise 
to perturbations that have the same gross observables as planetary perturbations. 

The severity of these degeneracies depends on both the specific par ameters of the plan- 
etary/ binar y companion , as we ll as on the data quality and coverage. IGriest fc Safizadeh 
(119981 ) and IChung et al.l (120051 ) demonstrated that the d <-> d~ l degeneracy is less severe 
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for more massive planets with separations closer to the Einstein ring (d ~ 1). Empiri- 
cally, this degeneracy was broken at the A% 2 ~ 4 level for the relatively large mass-ratio 
planetary companion OGLE-2005-BLG-071Lb (IDong et al.ll2008l ). for which the light curve 
was well-sampled, but was e ssentially unresolved for the low mass-ratio planetary compan- 
ion MOA-2007-BLG-192Lb (IBennett et al.ll2008l ). for which the pl anetary perturbation was 
poorly sampled. For the planetary /equal-mass binary degeneracy, lHan fc Gaudil (120081 ) ar- 
gued that, although the gross features of central caustic planetary perturbations can be 
reproduced by very close or very wide binary lenses, the morphologies differ in detail, and 
thus this degeneracy can be resolved with reasonable light curve coverage and moderate 



photometric precision. Indeed for every well-samp 
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the relatively poorly-sampled lig ht curve of MOA-20 07-BLG-192Lb, an equal-mass binary 
model is ruled out at A% 2 ~ 120 (IBennett et al.ll2008l ). although in this case this is partially 
attributable to the exquisite photometric precision (< 1%). 

One complication with searching for planets in high-magnification events is that, the 
higher the magnification, the more likely it is that the primary lens will transit the source. 
When this happens, the peak of the event is suppressed and smoothed out, as the lens 
strongly magnifies only a small portion of the source. If the source is also larger than the 
region of significant perturbations due to a planetary companion (roughly the size of the 
central caustic), then the planetary deviations will also be smoothed out and suppressed 
(jGriest Safizadehlll998l ; lHanll2007l ). These finite source effects have potential implications 
for both the detectability of central caustic perturbations, as well as the ability to uniquely 
determine the planetary parameters, and in particular resolve the two degeneracies discussed 
above. In practice, the caustic structures of all four high-magnification planetary events 
(containing five planets) were larger than the source. Hence, while there were detectable 
finite-source effects in all cases (which helped constrain the angular Einstein radius and 
so the physical lens parameters), the planetary perturbations were in all four cases quite 
noticeable. Thus the effect of large sources on the the detectability and interpretation of 
central-caustic perturbations has not been explored in practice. 

Theoretical studies of detectability of central caustic perturbati o ns when considerin g 
finit e sou r ce eff ects have been performed by lGriest fc Safizadehl (119981 ). IChung et al.l (120051 ). 
and iHanl (120071 ). These authors demonstrated that the qualitative nature of planetary per- 
turbations from central caustics is dramatically different for sources that are larger than 
the caustic. In particular, the detailed structure of the point-source magnification pattern, 
which generally follows the shape of the caustic, is essentially erased or washed out. Rather, 
the perturbation structure is characterized by a roughly circular region of very low-level, al- 
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most imperceptible deviation from the single-lens form that is roughly the size of the source 
and centered on the primary lens. This region is surrounded by an a nnular rim of larger 



deviations that has a w idth roughly equal to the width of the caustic fjGriest fc Safizadeh 



19981 ; IChung et al.ll2005l ). Finally, there are less pronounced deviations that extend to a few 
source radii. Planets are detectable even if their central caustics are quite a bit smaller than 
the source, provided that the \ 2 deviation is sufficiently high. (Often A% 2 > 60 is adopted, 
although Ax 2 > 150 may be more realistic.) The magnitude of these perturbations decrease 



(Chung et al. 


2005; 
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2007) 



Although central caustics may formally be detectable when the source is substantially 
larger than the caustic, it remains a significant question whether these very washed out 
caustics can be recognized in practice, and even if they can, whether they can be uniquely 
interpreted in terms of planetary parameters. Indeed, it is unknown whether a washed 
out central caustic due to a planet can actually be distinguished from one due to a binary 
companion. This question is especially important with regard to low-mass planets. The 
size of the central caustic scales as the product of the planet /star mass ratio and a definite 
function of planet-star separation. Hence, taken as a whole, smaller planets produce smaller 
caustics, meaning that events of higher magnification are required to detect them. These are 
just the events that are most likely to have their peaks washed out by finite-source effects. 

Here we analyze the first high-magnification event with a buried signature of a planet, 
in which the source size is larger than the central caustic of the planet. The caustic is indeed 
so washed out that the event appears unperturbed upon casual inspection. However, the 
residuals to a point-lens fit are clear and highly significant. We show that one can infer 
the planetary (as opposed to binary) nature of the perturbation from the general pattern 
of these residuals, and that a detailed analysis constrains the mass ratio of the planet quite 
well, but leaves the close/wide (d d~ l ) degeneracy intact. Hence, at least in this case, 
the fact that the caustic is buried in the source does not significantly hinder one's ability to 
uncover the planet and measure its mass ratio. 



2. Observations 

MOA-2007-BLG-400 [(a, <5) J20 oo.o = (18 h 09 m 41?98, -29°13'26.95"), (l,b)=(2.38 -4.70)] 
was announced as a probable microlensing event by the Microlensing Observations in Astro- 
physics (MOA) collaboration on 5 Sept 2007 (HJD' = HJD - 2450000 = 4349.1), about 5 days 
before peak. The source star proves to be a bulge subgiant and so is somewhat brighter than 
average, but the event timescale was relatively short (£e ~ 15 days) and observations had 
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been interrupted for 6 days by bad weather. Taken together, the two latter facts account for 
the relatively late alert. By coincidence, the triggering observations took place on the same 
night that another event, OGLE-2007-BLG-349 (aka MOA-2007-BLG-379), was peaking at 
extremely high magnification with an already-obvious planetary anomaly. After focusing 
exclusively on the latter event for the first 5 hours of the night, MOA resumed its normal 
field rotation for the last 1.5 hours, which led to the discovery of MOA-2007-BLG-400. 

The Microlensing Follow Up Network (//FUN) began observing this as a possible high- 
magnification event on 7 Sept, but did not mobilize intensive observations until UT 08:55 10 
Sept, just 15 hours before peak, following a high-mag alert issued by MOA a few min- 
utes earlier. Even at that point, the predicted minimum peak magnification was only 
Anax > 90, which would have enabled only modest sensitivity to planets. Nevertheless, 
all stops were pulled and it was observed as intensively as possible from 7 observatories, 
//FUN CBA Perth (Australia) 0.25m unfiltered, //FUN Bronberg (South Africa) 0.35m unfil- 
tered, //FUN SMARTS (CTIO, Chile) 1.3m V, J, H, //FUN Campo Catino Austral (CAO, 
Chile) 0.50m unfiltered, //FUN Farm Cove (New Zealand) 0.35m unfiltered, //FUN Auckland 
(New Zealand) 0.4m R, and //FUN Southern Stars (Tahiti) 0.28m unfiltered. 

The source star lies just outside one of the OGLE fields (as defined by their field tem- 
plates) and for this reason was not recognized as a microlensing event by the OGLE Early 
Warning System. However, due to small variations in pointing, there are a total of 452 
OGLE images containing this source. Only two of these are significantly magnified, ten days 
and nine days before peak. Hence, the OGLE data do not help constrain the light curve 
parameters. However, they are useful to study of the baseline behavior of the source (see 
Appendix [AT) . 

Essentially all of the "action" , both the peak of the event and the planetary anomalies, 
occurred during the //FUN SMARTS (Chile) observations at CTIO, 4354.47 < HJD' < 
4354.69, using the ANDICAM optical/IR dual-channel camera, and //FUN CAO (Chile) 
observations 4354.50 < HJD' < 4354.70. Most (45) of the optical CTIO observations over 
the peak were carried out in i band, with a few (8) taken in V in order to measure the (V — I) 
color. Each of these was a 5 minute exposure, with approximately 1 minute read-out time 
between exposures. During each optical exposure, there were 5 dithered if -band exposures, 
each of 50 seconds, almost equally spaced over the 6 minute cycle time. That is, 53 x 5 = 265 
if -band observations in all. Unfortunately, the source became so bright as it transited the 
lens (i.e., when the planetary anomalies were the strongest), that 14 i-band images were 
affected by non-linearities and saturation in the detector response. We exclude these 14 I- 
band data points from analysis. The if-band photometry are not affected by this problem, 
therefore, with higher time resolution and more continuous coverage than the i-band data, 
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the if -band data provide most of the constraining power to the microlens model. 

There were 84 CAO un- filtered observations taken during the peak night. Unfortunately, 
the peak of the event was severely saturated and the clock zero point is not securely known, 
therefore, these data are not used in the analysis. However, the exposures during the times 
of maximum deviation from a point lens (i.e., when the caustic was crossing the stellar limb) 
are not saturated, and these qualitatively confirm the interpretation from the more detailed 
CTIO data. 

The M OA data were reduced using the standard MOA differe nce imaging anal ysis (DIA) 



pipeline. All /iFUN data were reduced using DIA developed by I Wozniak! (120001 ). The H- 
band data are affected by intrapixel sensitivity variations at the 1% level. Fortunately, 
the dither pattern was repeated almost exactly over the night of the peak, so that these 
variations follow the 5-element dither pattern quite well. We therefore treat the if-band 
data as 5 independent data sets, which reduces y 2 by 180 for 8 degrees of freedom. The 
if-band images show a triangular PSF, which is likely to introduce systematic errors into 



the ph otometry. As a crosscheck, we also use the DIA package developed by iBond et al. 



(120011 ) to independently reduce these images. 



3. Microlens Model 

Despite the fact that the peak of MOA-2007-BLG-400 was "flattened" by finite-source 
effects, it nevertheless reached a very high peak magnification, A max = 628. However, even 
to the experienced eye, it looks like an ordinary point-lens light curve with pronounced finite- 
source effects. More detailed modeling is required to infer that it actually contains a Jovian 
mass-ratio planet. 

Figure [1] shows the light curve together with the best-fit point-lens model (blue) and 
planetary model (red). Both include finite-source effects. The most pronounced features 
of the point-lens model residuals are a short positive spike on the rising side and a short 
negative spike on the falling side, each lasting about 30 minutes, which leave very similar 
traces in I and H. As displayed in Figure [2], these features clearly stand out in the reductions 
using both the Wozniak (top panel) and Bond (bottom panel) DIA packages. Each package 
introduces its own systematic deviations, but there are no obvious trends besides the above 
features that are supported by both reductions. These occur very close to the times that the 
point lens begins and ends its transit of the source (within the framework of this model). The 
timing of these deviations strongly suggests that they are due to microlensing rather than 
stellar variability. There are then two possible explanations: either the source is actually 
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being transited by a more complicated caustic than a point lens (due to a binary or planetary 
companion) or the limb of the source is not being properly modeled. However, if one assumes 
a circular source, the latter explanation would imply symmetric residuals, whereas the actual 
residuals are closer to being antisymmetric. (We address the possibility of an elliptical source 
in Appendix |A]) Indeed, this approximate antisymmetry extends to the less pronounced 
residual features, including the sustained deficit prior to the first spike (and sustained excess 
following the second one) as well as the declining residuals between the two spikes. The 
similarity of the / and H residuals in itself argues that the deviations are due to microlensing 
rather than some sort of stellar variability, which would not generally be expected to be 
achromatic. 

The short durations of the spikes tell us that the central caustic is quite small, with 
"caustic width" w ~ 30min/15day = 1.4 x 10~ 3 . This implies that the companion either 
has low mass rat io, or is a very wide or very close binary companion. Formally, w is given by 



equation (12) of IChung et al.l (120051 ) as a very good approximation to the "short diameter" 
or "width" of a central caustic (see Fig. [3]). However, here the estimate is quite inexact not 
only because the width of the spikes is not precisely defined, but also because we do not 
know, at this point, the exact orientation of the caustic. 



After some algebra, one finds that in two limiting regimes, the IChung et al.l (120051 ) 
formula takes the forms, 

4(7 

w{d,q)^-£ (d»l), w(d,q) -> Aqd 2 (d < 1), (1) 

and 




where q is the companion/primary mass ratio and d is the separation in units of the Einstein 
ring. Note that in the first limit (eq. [I]), w — > 47, where 7 is the shear. The crossover point 
for these approximations is d — 2.3 (or d = 0.43), at which point each is in error by about 
15%. (For simplicity, we restrict the discussion here to the case d > 1. There is a well-known 
d <-> d^ 1 degeneracy between the (i<l and d 3> 1 limits, as can be guessed from the forms 
for w in eq. [T|. This degeneracy will prove to be almost perfect in this case, see § 13. 3 j) . Hence, 
in the first limit, d ~ 50g 1//2 , implying that if q were in the "binary range" (| log g| < 1), then 
d would be quite large. That is, the central caustic would be gener ated by a nearly pure 



shear and therefore would have a nearly symmetric, diamond-shaped, IChang-Refsdall (119791 ) 
form. In the point-lens model, the lens passes almost directly over the center of the source. 
For this trajectory, a symmetric caustic would yield symmetric residuals, in sharp contrast 
to Figure [H On the other hand, for the opposite limit, q ~ 2 x 10~ 3 (d — 1), which lies 
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squarely in the planetary regime. Thus, simple arguments already argue strongly in favor of 
a planetary companion that is fairly near the Einstein ring. 



3.1. Hybrid Pixel/Ray Map Algorithm 

Notwithstanding these arguments, we conduct a massive blind search for companions 
over a ve ry broad range of m asses using a modified version of the "magnification map" tech- 



nique of iDong et al.l (120061 ). which was specifically designed for high- magnification events. 
The original approach was, for each given (d, q) pair, to shoot rays over a fairly narrow an- 
nulus (say, 0.01 Einstein radii) around the Einstein ring in the image plane and to sort these 
rays in pixels on the source plane. Then for each source position being modelled (i.e., each 
data point), one would identify the pixels that intersected the source and would check each 
ray contained in these pixels to determine whether it landed on the source and, if so, evaluate 
the source surface brightness at that position. In the initial broad search, three parameters 
(d, q, a) are held fixed on a grid of values, while the remaining parameters (t , u , t E , p, and 
possibly others) are varied to minimize x 2 at each grid point. Here, a is the angle of the 
source trajectory relative to the binary axis, to is the time of closest approach to the adopted 
center of the lens geometry (usually the center of mass), uq is the source- lens separation at 
this time in units of the Einstein radius, p is the source radius in the same units, and t E is 
the Einstein crossing time. This division is efficient because 1) (to, Uo, ^e, p) are usually ap- 
proximately known from the general structure of the light curve, so minimization over these 
parameters is straightforward once (d, q, a) are fixed; 2) (d, q) define the map, which natu- 
rally facilitates minimization of other parameters except for a, whose value is not usually 
even approximately obvious from the light curve. 

The new approach differs principally in that the pixels that are contained entirely within 
the source are now evaluated as a whole, i.e., by the total number of rays in that pixel. Pixels 
that cross the source boundary are still evaluated ray-by-ray, as previously. This primary 
change then leads to several other changes. First, for each pixel, we record not only the 
geometric center but also the centroid of the rays. The surface brightness is then evaluated 
at the latter position. Second, the pixels are made much smaller, to minimize both the 
number of rays that must be evaluated individually and the surface-brightness variations 
across the pixel (which are corrected only to first order by the ray-centroid scheme just 
mentioned). Typically, there are a few hundreds of pixels per source. Third, the pixels are 
hexagonal, since this is the most compact tiling possible, i.e., the closest tile shape to a circle. 
Fourth, the source positions out side the map region ar e evaluated using the hex adecapole 
and quadrupole approximation of lGould J2008h ( see also lPejcha fc Heyrovsky||2007l ). Finally, 
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we use Markov chain Monte Carlo (MCMC) for the \ 2 minimization. 

3.2. The (w,q) Grid of Lens Geometries 

The initial search for solutions is conducted over a rectilinear grid in (w, q) rather 
{d,q). Since the short diameter, or "caustic width", w, is a monotonic function of the star- 
planet separation d (at fixed planet/star mass ratio q), these formulations are in some sense 
equivalent. However, for many events (including the present one) the short diameter w can 
be estimated by simple examination of the data. In these cases, the search space is both 
more regular and easier to define in terms of the (w, q) grid. In particular, equation (j2j) 
shows that at fixed w, d moves very close to 1 for very low q. 

3.3. Best-Fit Model 

We consider short diameters w over the range —3.5 < logw < —2 and companion 
mass ratios —4 < logg < 0, focusing on the regime d > 1. We find that there is only one 
local minimum in this range. The range of allowed solutions is well localized around this 
minimum, with 

q = 2.6 ±0.4 x 1(T 3 , d = 2.9 ± 0.2, w = 1.30 ± 0.06 x 10" 3 , (3) 

with the last quantity being, of course, dependent on the first two. Figure 0] shows the A% 2 = 
1, 4, 9 contours with respect to the mass ratio q and the projected planet-star separation d 
along with the short diameter w. The main point to note is that these parameters are quite 
well constrained. Note that, as expected, d and q are strongly correlated, while w and q 
are basically uncorrelated. We also perform similar searches using the alternative if-band 
reduction by Bond's DIA package. The solutions agree with the above to well within one 
sigma, but the parameters have larger uncertainties: q = 2.6 ± 0.7 x 10 -3 , d = 2.9 ± 0.3, w = 
1.26 ± 0.09 x 10 -3 . We therefore adopt results from the Wozniak-based reductions, noting 
that they may subject to systematic errors < la. 

Figure [3] shows the source trajectory and the central caustics as well as the differences in 
magnification between the best-fit planetary model and its corresponding single-lens model. 
This geometry nicely accounts for the main features of the point-lens residuals seen in Figure 
[TJ The regions beyond the "back walls" (long segments) of the caustic are somewhat de- 
magnified, which accounts for the initial depression of the light curve. As the source crosses 
the "back wall" of the caustic, it spikes. After the source has exited the caustic, it continues 
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to suffer additional magnification due to the "ridge" of magnification that extends from the 
trailing cusp. 

We also conducted a similar blind search as above, but concentrating on the regime 
d < 1. As expected, we recover the well-known d «-> d^ 1 degeneracy, and find a solution with 
essentially the same q, but with d = 2.9 -1 = 0.34 ± 0.02, and the wide solution is slightly 
preferred by A% 2 = 0.2. Thus, although each solution is well-localized to its respective 
minimum, this discrete degeneracy implies that the projected separation can take on two 
values that differ by a factor of ~ 8.5. The severity of the degeneracy can be traced to 
the planetary parameters. Although the planet/star mass ratio is quite large, which tends 
to reduce the severity of the degeneracy, the planet lies quite far from the Einstein ring, 
which tends to make it more severe. Actually, a better measure of the overall expected 
asymmetry between the d and d~ x solutions is the short diameter w, which in this case 
is small, implying a severe degeneracy. Indeed, the caustic structure and magnification 
pattern of the two solutions are nearly identical. In this case, the large size of the source has 
competing influences on the ability to resolve the degeneracy. On one hand, the large size of 
the source serves to suppress the planetary deviations, thus making subtle differences more 
difficult to distinguish. On the other hand, the large source implies that a large fraction 
of the planetary perturbation region is probed. In this case, the source probes essentially 
the entire region of significant planetary perturbation, as can be seen in Figure [31 This 
is important for distinguishing between the solutions, as the largest difference between the 
magnification patterns of th e two degenerate solutio ns occurs in the region near the tip 



of the arrow-shaped caustic (IGriest fc Safizadehlll998l ). From Figure [3] it is clear that this 



region would have been entirely missed if the source had been substantially smaller than the 
caustic. 



4. Finite-Source Effects 

In addition to (d, q), the model also yields the source radius relative to the Einstein 
radius, 

p = 9*/9 E = 3.29 ± 0.08 x 10" 3 , (4) 



We then follow the standard (jYoo et al.l 120041 ) technique to determine the angular source 
radius, 

0* = 1.05 ± 0.05 /xas. (5) 

That is, we first adopt [(V — I)o, Io]cium P = (1.00, 14.32) for the dereddened position of the 
clump. We then measure the offset of the source relative to the clump centroid A[(V—I), I] = 
(-0.19,3.25), to obtain [(V - I) ,io]« = (0.81,17.57). See Figured The instrumental 
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source color is derived from model-independent regression of the V and I flux, while the 
instrumental magnitude is obtained fro m the light-curve mode l. We convert (V — I) to 



(V — K) using the color-color relations of iBessell fc Brettl (119881 ) . yielding (V — K) n = 1.75, 



and th en obtain equation using the color/surface-brightness relations of iKervella et al. 



M = 0.10 ± 0.01 M ( ^ ) (6) 



( 120041 ) . Combining equations (@J and (jSJ) gives 6e = 0*/ p = 0.32 mas. And combining this 
with the definition 0| = K,Mii ve \, where M is the lens mass, 7r re i is the source-lens relative 
parallax, and k = 4G/c 2 AU ~ 8.1 masMg 1 , together with the measured Einstein timescale, 
£e = 14.3 ± 0.3days, we obtain 

r = 0.10 ± 0.01 mA 

. 125 yuas 
and 

/i re i = — = 8.2 ± 0.5masyr~ . (7) 

The relatively high lens-source relative proper motion /i re i is mildly suggestive of a foreground 
disk lens, but still quite consistent with a bulge lens. Since 7r re i = 125 /ias corresponds to a 
lens distance D L = 4kpc (assuming source distance D s = 8kpc), equation (jSD implies that 
if the lens did lie in the foreground, then it would be a very low-mass star or a brown dwarf. 

Assuming that the source lies at a Galactocentric distance modulus 14.52, its dereddened 
color and magnitude imply that [(V — I)q, Mj] = (0.81, 3.07), making it a subgiant. 



5. Limb Darkening 



As illustrated in Figure [H the principal deviations from a point-lens light curve occur 
at the limb of the star. This prompts us to investigate the degree to which the planetary 
solution is influenced by our treatment of limb darkening. The results that we report are 
based on a fit to the if -band surface brightness profile of the form 



m 

On 



r i 



3 9 

- cosv 



Al 



^cos 1 / 2 ^ 



where T and A are the linear and square-root parameters, respectively, and where d is the 
angle between t he normal to the stellar surface and the line of sight, i.e., sin^ = See 
An et al.l (120021 ) for the relation between (r, A) and the usual (c, d) formalism. 



In deriving the reported results, we fix the if -band limb- darken i ng par ameters (r, A) = 
(—0.15,0.69), corresponding to (c,d) = (—0.21,0.79) from Oaretl (120001 ) for a star with 
effective temperature T e s = 5325 K and lo g q = 4.0. These stellar parameters are suggested 
by comparison to Yale-Yonsei isochrones (jDemarque et al.l 120041 ) for the dereddened color 
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and absolute magnitude reported in § HI We also perform fits in which F and A are allowed to 
be completely free. From these fits, we find that our best-fit model has (T, A) = (—0.64, 1.47). 
T and A are highly correlated, so their individual values are not of interest, and the surface- 
brightness profiles generated by these two sets of (r, A) are qualitatively similar. In the 
presen t context, how ever, the key point is that when we fix the limb-darkening parameters 
at the IClaretl (120001 ) values, the contours in Figure |4] remain essentially identical and the 
best fit values change by much less than 1 a. 

Because of lower point-density and the aforementioned problems with the / data over 
the peak, we only attempt a linear l imb-d arkening fit, i.e., we use equation ([8]) with A = 0, 
and we adopt V = 0.47 from|ciaret| (hoooh . 



6. Blended Light 

In the crowded fields of the Galactic bulge, the photometered light of a microlensing 
event rarely comes solely from the lensed source. Rather there is typically additional light 
that is blended with the source but is not being lensed. This light can arise from unrelated 
stars that happen to be projected close enough to the line of sight to be blended with the 
source, or it can come from companions to the source, companions to the lens, or the lens 
itself. This last possibility is most interesting because, if the lens flux can be isolated and 
measured, it provides strong constraints on the lens properties, and in this case would enable 
a compl ete solution of the le ns mass and distance, when combined with the measurement of 



6 E ( e.g. . Bennett et all 120071 ). 



To investigate the blended light, we begin by using the method of iGould fc Anl (120021 ) 
to construct an image of the field with the source (but not the blended light) removed, 
and compare this to a baseline image, which of course contains both the source and the 
blended light. In these images, the source/blend is immersed in the wings of a bright star 
(roughly 3.7 mag brighter than the source), which lies about 2" away. On the baseline 
image, the source/blend is noticeable against this background, but hardly distinct. On the 
source-subtracted image, the blend is not directly discernible. 

To make a quantitative estimate of the blend flux, we fit the region in the immediate 
vicinity of the bright star to the form F = ai + a 2 x PSF, where "PSF" is the point-spread- 
function determined from the DIA analysis. We then subtract the best fit flux profile from 
the image. This leaves a clear residual at the position of the source/blend in the baseline 
image, but just noise in the source-subtracted image. We add all the flux in a 1.8" square 
centered on the lens, finding 564 ADU and —28 ADU, respectively. We conclude that the 
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J-band flux blend/source ratio is fbj f s < 0.05. Of course, even if we had detected blended 
light, it would be impossible to tell whether it was directly coincident with the source. If it 
were, this would imply that this light would be directly associated with the event, i.e., being 
either the lens itself or a companion to the lens or the source. Hence, this measurement is an 
upper limit on the light from the lens in two senses: no light is definitively measured, and if it 
were we do not know that it came from the lens. Combining this limit with equation (JBJ), and 
assuming the lens is a main-sequence star, it must then be less massive than M < 0.75 M , 
and so must have relative parallax 7r re i > 15 mas. This implies a lens-source separation 
Ds — Dl > 1 kpc, which certainly does not exclude bulge lenses. Indeed, if the lens were a 
K dwarf in the Galactic bulge, it would saturate this limit. 



7. Discussion 

MOA-2007-BLG-400 is the first high-magnification microlensing event for which the 
central caustic generated by a planetary companion to the lens is completel y enveloped by 



the s ource. As a comparison, the planetary caustic of OGLE-2005-BLG-390 (IBeaulieu et al. 



20061 ) is smaller than its clump-giant source star in angular size. When the planetary caustics 



is covered by the source, the finit e-source effects broade n the "classic" iGould fc Loebl (Il992l ) 



planetary perturbation features (IGaudi &: Gouldlll997l ). By contrast, planet-induced devia- 



tions in MOA-2007-BLG-400 are mostly obliterated, rather than being broadened, because 
the source crosses the central caustic rather than the planetary caustic. We showed, never- 
theless, that the planetary character of the event can be inferred directly from the light-curve 
features and that the standard microlensing planetary parameters (d, q) = (2.9, 2.6 x 10~ 3 ) 
can be measured with good precision, up to the standard close/wide d <-> d^ 1 degeneracy. 
We demonstrated that, in this case, the close/wide degeneracy is quite severe, and the wide 
solution is only preferred by A% 2 = 0.2. This is unfortunate, since the separations of the 
two solutions differ by a factor of ~ 8.5. We argued that the severity of this degeneracy 
was primarily related to the intrinsic parameters of the planet, rather than being primarily 
a result of the large source size. 

Although the mass ratio alone is of considerable interest for planet formation theories, 
one would also like to be able to translate the standard microlensing parameters to physical 
parameters, i.e., the planet mass m p = qM, and planet-star projected separation r±_ = 
dO-ftD^. Clearly this requires measuring the lens mass M and distance Dl. In this case, 
the pronounced finite source effects have already permitted a measurement of the Einstein 
radius #e = 0.32 mas, which gives a relation between the mass and lens-source relative 
parallax (eq. [6]). This essentially yields a relation between the lens mass and distance, since 



-17- 



the source distance is close enough to the Galactic center that knowing Dl is equivalent to 
knowing 7r re i. Therefore, a complete solution could be determined by measuring either M or 
Dl, or some combination of the two. 

One way to obtain an independent relation between the lens mass and distance is to 
measure the microlens parallax, tte- There are two potential ways of measuring 7Te- First, 
one can measure distortions in the light curve arising from the acceleration of the Earth as it 
moves along its orbit. Unfortunately, this is out of the question in this case because the event 
is so short that these distortions are immeasurably small. Second, one can measure the effects 
of terrestrial parallax, which gives rise to differences between the light curves simultaneously 
observed from two or more observatories separated by a significant fraction of the diameter 
of the Earth. Practically, measuring these differences requires a high-magnification event, 
which would appear to make this event quite promising. Unfortunately, although we obtained 
simultaneous observations from two observatories separated by several hundred kilometers 
during the peak of the event, one of these datasets suffers from large systematic errors and 
an unknown time zero point, rendering it unusable for this purpose. 

The only available alternative for breaking the degeneracy between the lens mass and 
distance would be to measure the lens flux, either under the glare of the source or, at a later 
date, to separately resolve it after it has moved away from the line of sight to the source 



( lAlcock et al.ll200ll ; iKozlowski et al.ll2007l ). Panels (e) and (f) of Figure [6] show the Bayesian 
estimates of the lens brightness in /-band and if -band, respectively. If the lens flux is at least 
2% of the source flux, then the former kind of measurement could be obtained from a single 
epoch Hubble Space Telescope observation, provided it were carried out in the reasonably 
near future. At roughly 99% probability, the blended light would be either perfectly aligned 
with the source (and so associated with the event) or well separated from it. HST images 
can be photometrically aligned to the ground-based images using comparison stars with 
an accuracy of better than 1%. Hence, photometry of the source+blend would detect the 
blend, unless it were at least 4 mag fainter than the source. In principle, the blend could be 
a companion to either the source or lens. Various argumen ts can be us e d to c onstrain either 



of those scenarios. We do not explore those here, but see iDong et all (120081 ). If the lens is 
not detectable by current epoch HST observations (or no HST observations are taken), then 
it will be detectable by ground-based AO if -band observations in about 5 years. This is 
because the lens-source relative proper motion is measured to be /i re i = 8masyr _1 , and the 
diffraction limit at H band on a 10m telescope is roughly 35 mas. If the lens proves to be 
extremely faint, then a wider separation (and hence a few years more time baseline) would 
be required. 

In the absence of additional observational constraints, we must rely on a Bayesian anal- 
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ysis to estimate the properties of the host star and p lanet, which incorporates priors on the 



distribution of lens masses, distances, and velocities (lDominikll2006l ; iDong et al.1120061). This 



i s a s tandard procedure, which we only briefly summarize here. We adopt a lHan fc Gould 



( 119951 ) model for the Galactic bar, a double-exponential disk with a scale height of 325 pc, 



and a scale leng t h of 3.5 kpc, as well as other Galactic model parameters as described in 



Bennett fc Rhid (120021 ). We incorporate constraints from our measurement of the lens an- 
gular Einstein radius #e and the event timescale, as well as limits on the microlens parallax 
and /-band magnitude of the lens. In practice, only the measurements of 9e and provide 
interesting constraints on these distributions. In addition, we include the small penalty on 
the close solution, exp(— Ax 2 /2), where the wide solution is favored by A% 2 = 0.2. For the 
estimates of the planet semimajor axis, we assume circular orbits and that the orbital phases 
and cos(inclinations) are randomly distributed. 

The resulting probability densities for the physical properties of the host star, as well 
as selected properties of the planet, are shown in Figure El The Bayesian analysis suggests 
a host star of mass M = O.3Ol o 'i2^0 & t distance of Dl = 5. 8^8 kpc. In other words, 
given the available constraints, the host is most likely an M-dwarf, probably in the fore- 
ground Galactic bulge. Given that the planet/star mass ratio is measured quite precisely, 
the probability distribution for the planet mass is essentially just a rescaled version of the 
probability distribution for the host star mass. We find m p = 0.82^33 Mj up . The close/wide 
degeneracy is apparent in the probability distribution for the semimajor axis a. We estimate 
^ciose = 0.72l°;fg AU for the close solution, and a w id e = 6.5 jlf ;| AU for the wide solution. The 
equilibrium temperatures for these orbits are T eq close = 103^26 K anc ^ ^eq.,wide = 34 ± 9 K 
for the close and wide solutions, respectively. 

Thus our Bayesian analysis suggests that this system is mostly likely a bulge mid- 
M-dwarf, with a Jovian-mass planetary companion. The semimajor axis of the planetary 
companion is poorly constrained primarily because of the close/ wide degeneracy, but the 
implied equilibrium temperatures are cooler than the condensation temperature of water. 
Specifically we find that T eq < 173 K at 2 a level. Alternatively, if we assume the snow 
line is given by a snow = 2.7 AU(M/M ), we find for this system a snow line distance of 
~ 0.81 AU, very close to the inferred semimajor axis of the close solution. Thus this planet 
is quite likely to be located close to or beyond the snow line of the system. 

Although we cannot distinguish between the close and wide solutions for the planet 
separation, theoretical prejudice in the context of the core-accretion scenario would suggest 
that a gas-giant planet would be more likely to form just outside the snow line, thus preferring 
the close solution. However, we have essentially no observational constraints on the frequency 
and distribution of Jupiter-mass planets at the separations implied by the wide solution 
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(~ 5.3 — 9.7 AU), for such low-mass primaries. Unfortunately, the prospects for empirically 
resolving the close/wide degeneracy in the future are poor. The only possible method of doing 
this would be to measure the radial velocity signature of the planet. Given the faintness of 
the host star (see £0 and Fig. [6]), this will likely be impossible with current or near-future 
technology. 

The mere existence of a gas-giant planet orbiting a mid-M-dwarf is largely unexpected 
in the core-accretio n scenario, as formatio n of such planets is thought to be inhibited in such 
low-mass primaries (ILaughlin et al.ll2004l ). Observationally, however, although the frequency 
of Jovian companions to M-dwarfs with a < 3 AU does appear to be smaller than the 



corre s ponding frequency of s uch companions to FGK dwarfs (jEndl et al.ll2006l ; iJohnson et al. 



2007 



Cumming et al.ll2008l ). several Jovian- mass companions to M dwarfs are known (see 



Dong et al.ll2008l for a discussion), so this system would not be unprecedented. Furthermore, 
it must be kept in mind that the estimates of stellar (and so planet) mass depend on the 
validity of the priors, and even in this context have considerable uncertainties. 

Most of the ambiguities in the interpretation of this event would be removed with a 
measurement of the host star mass and distance, which could be obtained by combining our 
measurement of #e with a measurement of the lens light as outlined above. The Bayesian 
analysis informs the likelihood of success of such an endeavor. This analysis suggests that, 
if the host is a main-sequence star, its magnitude will be II = 23.9l^o anc ^ Hl = 21.4^0, 
which corresponds to 0.6% and 1.7% of the source flux, respectively. If initial efforts to detect 
the lens fail, more aggressive observations would certainly be warranted: microlensing is the 
most sensitive method for detecting planets around very low-mass stars simply because it is 
the only method that does not rely on light from the host (or the planet itself) to detect the 
planet. And given equation ([6]), even an M dwarf at the very bottom of the main sequence 
M = 0.08 M Q , would lie at D L = 3.5 kpc and so would be H ^ 24. 
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A. Failure of Elliptical-Source Models 

Because MOA-2007-BLG-400 is the first microlensing event with a completely buried 
central caustic, it is important to rule out other potential causes of the deviations seen in 
the light curve (apart from a planetary companion to the lens). The principal features of 
these deviations are the twin "spikes" in the residuals, which are approximately centered on 
the times when the lens enters and exits the source. In the model, these crossings occur at 
about HJD' = 4354.53 and HJD' = 4354.63, i.e., very close to the spikes in Figured! In 
principle, one might be able to induce such spikes by displacing the model source crossing 
times from the true times. The only real way to achieve this (while still optimizing the 
overall fit parameters) would be if the source were actually elliptical, but were modeled as a 
circle (which, of course, is the norm). 

One argument against this hypothesis is the similarity of the / and H residuals (§[3]). If 
the source were an ellipsoidal variable, then one would expect color gradients due to "gravity 
darkening" . 

Nevertheless, we carried out two types of investigation of this possibility. First, we 
modeled the light curve as an elliptical source magnified by a point (non-binary) lens. In 
addition to the linear flux parameters (source flux plus blended flux for each observatory) 
there are 6 model parameters, the three standard point-lens parameters (to,uo,t^,), plus the 
source semimajor and semiminor axes (p a , pb) and the angle of the source trajectory relative 
to the source major axis, a. We find that the elliptical source reduces x 2 by about 200, 
but it does not remove the "spikes" from the residuals, which was the primary motivation 
for introducing it. Instead, essentially all of the \ 2 improvement comes from eliminating 
the asymmetries from the rest of the light curve. Recall, however, that the planetary model 
removes both these asymmetries and the "spikes" . Moreover, the best-fit axis ratio is quite 
extreme, Pb/p a = 0.7, which would produce very noticeable ellipsoidal variations unless the 
binary were being viewed pole on. 

Next we looked for sinusoidal variations in the baseline light curve. The individual 
OGLE errorbars at baseline are smaller than for MOA, and since ellipsoidal variations are 
strictly periodic, the longer OGLE baseline (about T = 2000 days versus T = 800 days for 
MOA) does a better job of isolating this signal from various possible systematics. Therefore, 
for this purpose, the OGLE data are more suitable than MOA. The OGLE data are essentially 
all baseline (only two magnified points out of 452). Their periodogram shows several spikes 
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at the 0.01 mag level, and a maximum A% 2 = 20, which are consistent with noise. The width 
of the spikes is extremely narrow, consistent with the theoretical expectation for uniformly 
sampled data of a(P)/P 2 ~ a/24/Ax 2 /(27tT) ~ 10 _4 day _1 , indicating that the data set is 
behaving normally. 

In brief, our investigation finds no convincing evidence for ellipticity of the source, 
certainly not for the several tens of percent deviation from circular that would be needed 
to significantly ameliorate the deviations seen near peak in the light curve. Moreover, even 
arbitrary source ellipticities cannot reproduce the light curve's most striking features: the 
two "spikes" in the residuals that occur when the lens crosses the source boundary. 
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Fig. 1.— Top: Lightcurve of MOA-2007-BLG-400 with data from /xFUN CTIO (Chile) si- 
multaneously taken in H (cyan), I (DIA, black). Models are shown for a point lens (blue) 
and planet-star system (red). There are 5 50-second H exposures for each 300-second / (or V 
- not shown) exposure in 6 minutes cycles. Some /-band data at the peak suffer from satura- 
tion, and those points are therefore removed from the analysis (see text). Middle: Residuals 
for best-fit point-lens model and its difference with the planetary model. Note that in the 
top panel, the H data are shown as observed, while the / data are aligned. Normally, such 
alignment is straightforward because microlensing of point sources is achromatic. However, 
here there is significant chromaticity due to different limb-darkening. The /-band points in 
the top panel are actually the residuals to the /-band limb-darkened model (middle panel), 
added to the iZ-band model curve (top panel). Bottom: Residuals from a point-lens model 
with the same parameters as the planetary model, which can be directly compared to the 
"magnification map" in Fig. [3j These "didactic residuals" are naturally more pronounced 
than those from the best-fit point lens. 
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Fig. 2. — Comparisons of residuals to the best-fit point-lens models between two photometric 
reduction s of /xFUN i f -ban d data, using the DIA packages developed by IWozniak! (120001 ) 
(top) and lBond et al.l (120011 ) (bottom). The red curve in each panel represents the difference 
of the best-fit planetary and point-lens models for that panel's reduction. Both reductions 
agree on the main planetary features, but each package introduces its own systematics. 
For example, the systematic deviations from the planetary model at HJD ~ 2454354.57 
shown in the top panel are not supported by the reductions of the Bond's package. During 
2454354.46 < HJD < 2454354.51, most images have low transparency (< 50%), which 
causes relatively large scatter in Bond's DIA reductions. In comparison, the reduction by 
Wozniak's DIA has smaller scatter during this period. However, the data exhibit some 
low-level systematics, which are not supported by the other reduction. 
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Fig. 3. — Magnification differences between of best-fit planetary model [(q, d) = (0.0026, 2.9) 
and (q, d) = (0.0026, 0.34) being nearly identical] and single-lens models, in units of the 
measured source size, p = 0.0033 Einstein radii. Contours show 1%, 2%, 3%, and 4%, 
deviations in the positive (brown) and negative (blue) directions. Top panel: Single-lens 
geometry (t , u , £g) is taken to be the same as in the planetary model, with no finite source 
effects. Caustic (contour of infinite magnification) is shown in white. The deviations are very 
pronounced. Bottom Panel: Same as top panel, but including finite-source effects, which 
now explain the main features of the light curve. The trajectory begins with a negative 
deviation, then hits a narrow "brown ridge" causing the spike seen in the bottom panel of 
Fig. [H as the edge of the source first hits the caustic. Then there are essentially no deviations 
(white) while the source covers the caustic. The caustic exit induces a narrow "blue ridge" 
corresponding to the negative-deviation spike seen in Fig [JJ Finally, the source runs along 
the long "brown ridge" corresponding to the prolonged post-peak mild excess seen in Fig [TJ 
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Fig. 4. — Contours of A% 2 = 1, 4, 9 relative to the minimum as a function of planet-star mass 
ratio q and projected planet-star separation d (top), as well as "short caustic diameter" (see 
Fig. [3]) w (bottom), w (in units of 9e) is a function of q and d (see text). The solution shown 
here corresponds to q = 2.6 ± 0.4 x KT 3 and d = 2.9 ±0.2 (or d = 0.34 ±0.02). These values 
of d correspond to physical separations and equilibrium temperatures of ~ 0.6 — 1.1 AU, 
~ 103 K and ~ 5.3 — 9.7 AU, ~ 34 K for the close and wide solutions, respectively. 



(V-l) [instrumental] 



Fig. 5. — Instrumental color-magnitude diagram of field containing MOA-2007-BLG-400. 
The color and magnitude of the source (black) are derived from the fit to the light curve, which 
also yields an upper limit for the J-band blended flux (green) . The large error bar on the latter 
point indicates a complete lack of information about its \^-band flux. The clump centroid 
is indicated in red. From the source-clump offset, we estimate [I, (V — J)]o, a = (17.57, 0.81), 
implying it has angular radius 6* = 1.05 //as. Assuming the source lies at 8 kpc, it has 
[Mj, (V — I)]o,« = (3.07, 0.81), making it a subgiant. The lack of blended light allows us to 
place an upper limit on the lens flux, which implies that it has mass M < 0.75 M & . See text. 
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M0A-2007-BLG-400Lb Properties 
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Fig. 6. — Bayesian relative probability densities for the physical properties of the planet 
MOA-2007-BLG-400Lb and its host star, (a) Mass of the host star, (b) Planet semimajor 
axis, (c) Distance to the planet/star system, (d) Equilibrium temperature of the planet, (e) 
/-band magnitude of the host star, (f) if -band magnitude of the host star. In panel (a), 
we also show the probability density for the planet mass, which is essentially a rescaling of 
that of the host star, because the mass ratio is measured so precisely q = 2.6 ± 0.4 x 10~ 3 . 
In all panels, the solid vertical lines show the medians, and the 68.3% and 95.4% confidence 
intervals are enclosed in the dark and light shaded regions, respectively. In panel (b) and (d), 
the probability distributions for wide and close degenerate solutions are computed separately, 
and then the wide solution is weighted by exp(— A% 2 /2), where A% 2 = 0.2 is the difference 
between them. These distributions are derived assuming priors obtained from standard 
models of the mass, velocity, and density distributions of stars in the Galactic bulge and 
disk, and include constraints from the measurements of lens angular Einstein radius 6*e 
and the timescale of the event £e, as well as limits on the /-band magnitude of the lens. 
In practice, only the measurements of #e and £e provide interesting constraints on these 
distributions. 
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Table 1. Best-fit Planetary Models 



Model 1 


to — t m { 2 


w 


*E 


d 


Q 


a 3 


P 




day 




day 






deg 




Close 


0.08107 


0.00025 


14.41 


0.34 


0.0026 


227.06 


0.00326 


Wide 


0.08106 


0.00027 


14.33 


2.87 


0.0025 


226.99 


0.00329 



1 The wide solution is preferred over best-fit single-lens model by A% 2 = 
1070.04 and the close solution by Ax 2 = 1069.84 

2 t rcf = HJD 2454354.5 

3 The geometry of the source trajectory is illustrated in Figure [31 in 
which the planet is to the right of the lens star, to, uq, and a are defined 
with respect to the "center of magnification", which is the center of mass 
of the star/planet system for the close model and q/(l + q)/d away from 
the position of the lens star toward the direction of the planet for the wide 
model. 



